Quantum state detection of a superconducting flux qubit using a DC-SQUID in the 

inductive mode 
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We present a readout method for superconducting flux qubits. The qubit quantum flux state 
can be measured by determining the Josephson inductance of an inductively coupled DC super- 
conducting quantum interference device (DC-SQUID). We determine the response function of the 
DC-SQUID and its back-action on the qubit during measurement. Due to driving, the qubit energy 
relaxation rate depends on the spectral density of the measurement circuit noise at sum and differ- 
ence frequencies of the qubit Larmor frequency and SQUID driving frequency. The qubit dephasing 
rate is proportional to the spectral density of circuit noise at the SQUID driving frequency. These 
features of the backaction are qualitatively different from the case when the SQUID is used in the 
usual switching mode. For a particular type of readout circuit with feasible parameters we find that 
single shot readout of a superconducting flux qubit is possible. 

PACS numbers: 03.67.Lx, 03.65.Yz, 85.25.Cp , 85.25.Dq 



I. INTRODUCTION 

An information processor based on a quantum mechan- 
ical system can be used to solve certain problems signif- 
icantly faster than a classical computeri This idea has 
motivated intense research in recent years on the control 
and measurement of quantum mechanical systems. The 
basic units in a quantum computer are two level systems, 
also called quantum bits or qubits. Many types of qubits 
based on various physical systems have been proposed 
and implemented experimentally. 

Qubits based on solid state systems have the advan- 
tage of flexibility in design parameters and scalability. 
An important class of solid state qubits are the super- 
conducting qubits. They are mesoscopic systems formed 
of superconductor structures containing Josephson junc- 
tions. The energy level structure in these systems is the 
result of the interplay between the charging energy, as- 
sociated with the electrostatic energy due to distribu- 
tion of the charge of a single Cooper pair, and of the 
Josephson energy, associated with the tunnelling prob- 
ability for Cooper pairs across the Josephson junctions. 
Quantum coherent oscillations have been observed for a 
few versions of qubits with Josephson junctions2iiMi£i§i£ 
and coupling of two qubits was demonstrated—. 

A suitable qubit state detection apparatus for individ- 
ual qubits is an essential ingredient for the implemen- 
tation of algorithms for a quantum computer. Efficient 
measurement is necessary to extract all the relevant in- 
formation on single qubit states within a restricted time. 
Moreover, for correlation type measurements in a mul- 
tiple qubit system, the unwanted backaction of the first 
measurement should not disturb the system so strongly, 
that subsequent measurements will be meaningless. In 
this paper we discuss a measurement method for super- 
conducting flux qubits. Flux qubits are a qubit variety 
formed of a superconducting loop interrupted by Joseph- 
son junctions. The basis states have oppositely circulat- 



ing persistent currents in the loop. The control param- 
eter is an external magnetic flux in the qubit loop. The 
qubit state can be determined by measuring the magnetic 
flux generated by its persistent current. 

A natural candidate for the measurement of the state 
of a flux qubit is a DC superconducting quantum in- 
terference device (DC-SQUID). A DC-SQUID is a loop 
containing two Josephson junctions. Its critical current, 
which is the maximum supercurrent that it can sustain, 
depends on the magnetic flux enclosed in the loopS. The 
state of a flux qubi^Siii, was measured using an un- 
derdamped DC-SQUIDLiS. The critical current of the 
SQUID, and thus the state of the flux qubit, is deter- 
mined as the maximum value of the current, where the 
SQUID switches to a finite voltage state. Due to thermal 
and quantum fluctuations, switching of the SQUID is a 
stochastic process^. The qubit states are distinguish- 
able if the difference between the two average values of 
the switching current, corresponding to the qubit flux 
states, is larger than the statistical spread of the mea- 
sured values of the switching current. The measurement 
of a flux qubit using a switching DC-SQUID was char- 
acterized by an efficiency as large as 60 %Z. Further 
improvement of the measurement efficiency is possible 
and the backaction on the measured qubit if no measure- 
ment is performed can be reduced to acceptable levels^. 
Nevertheless, switching to the dissipative state has a few 
drawbacks. The finite voltage across the DC-SQUID de- 
termines the generation of quasiparticles which causes 
decoherence of the qubit^. The long quasiparticle re- 
combination time is a severe limit to the reset times for 
the qubits. In the finite voltage state the SQUID gener- 
ates AC signals with frequencies in the microwaves range 
and broad spectral content, that can induce transitions 
in a multiple qubit system, constrained to have energy 
levels spacings in the same region. The mentioned types 
of backaction will not have an affect on the statistics of 
the measurements on a single qubit; however, in a com- 
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plex multiple qubit system switching of a DC-SQUID to 
the finite voltage state is a strong disturbance of the state 
of the total system which introduces errors in subsequent 
computations and/or measurements. 

A DC-SQUID can be used as a flux detector in an 
alternative mode of operation, in which switching to the 
dissipative state is avoided. This is based on the property 
of a SQUID to behave as an inductor, with a Josephson 
inductance that depends on the magnetic flux enclosed in 
the loorA The value of the Josephson inductance can be 
determined by measuring the impedance of the SQUID. 
The flux sensitivity in this operation mode is increased 
if the junction is shunted by a capacitor and the circuit 
is excited with an AC signal at a frequency close to the 
resonance frequency. A SQUID in the inductive mode 
integrated in a resonant circuit was used for the mea- 
surement of spectroscopy of a flux qubifei^. 

The inductive operation mode resembles the RF- 
SQUID in the dispersive mode 9 . The RF-SQUID con- 
tains a superconducting loop with a single Josephson 
junction; the impedance of a high quality tank circuit in- 
ductively coupled to the loop is measured near resonance, 
where it is very sensitive to the value of the magnetic flux 
in the loop. For charge measurement, a similar device is 
the RF single electron tunnelling transistor (RF-SET), 
with the difference that a dissipative property of a SET 
transistor is measured directly. The RF-SET was used 
as a detector for charge qubits by Duty et al 6 . Moti- 
vated by research on superconducting qubits a few flux 
or charge detectors based on the measurement of a reac- 
tive circuit element have been recently implemented. A 
flux qubit was studied by Grajcar et ali± by measuring 
the susceptibility of the qubit loop using a coupled high 
quality tank circuit. A detector for charge qubits based 
on the measurement of the inductance of a superconduct- 
ing SET was proposed by ZoriniSi and implemented by 
Sillanpaa, et alm&. A sensitive measurement of the critical 
current of a Josephson junction which exploits the non- 
linearity of the current phase relation was demonstrated 
by Siddiqi et alM. The state of a charge qubit was read 
out by Wallraff et alm^ by measuring the transmission 
through a coupled transmission line resonator. 

The paper is organized as follows. In Sec. [H]we dis- 
cuss a few general constraints on the parameter range 
where the DC-SQUID in the inductive mode can oper- 
ate. We continue in Sec. IIIII with a general analysis of 
the response function of this device as a flux detector. 
The response function is derived for a general type of 
circuit embedding the DC-SQUID. In Sec. EJwe discuss 
the qubit-SQUID interaction and we identify the relevant 
aspects of the measurement backaction. The energy re- 
laxation rate and the dephasing rate of the qubit during 
the measurement are derived in Sec. E] Because of AC 
driving of the SQUID and quadratic coupling of the qubit 
to the SQUID, the qubit relaxation rate is proportional to 
the spectral density of circuit noise at frequencies which 
are the sum and the difference of SQUID AC driving 
frequency and qubit Larmor frequency. Similarly, the 
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FIG. 1: (a) Ground and excited state energy levels for a PCQ, 
with Ej = 258 GHz, E c = 6.9 GHz and a = 0.75, which 
results in I p = 300 nA and A = 5.5 GHz. (b) Expectation 
value of the loop current for the ground and excited energy 
states with the parameters mentioned in (a), (c) Schematic 
representation of the PCQ and of the measuring DC-SQUID, 
with crosses indicating Josephson junctions. The SQUID acts 
like a variable inductor, with a value dependent on the qubit 
flux state. The impedance of a resonant circuit formed of the 
DC-SQUID and a shunt capacitor depends on the qubit state 
(R represents an equivalent series resistance for the resonance 
circuit, not shown). 



dephasing rate is proportional to the spectral density of 
circuit noise at the frequency of the SQUID AC driving. 
In Sec. IVII we discuss the results of the calculations on 
the measurement backaction. We analyze the measure- 
ment efficiency for a specific readout circuit and we find 
that single shot readout of a flux qubit is possible. 



II. GENERAL CONSIDERATIONS 

In this paper we focus on the readout of a persistent 
current qubit (PCQ)i2*ii, though the analysis of the mea- 
surement is applicable to flux qubits in general. The PCQ 
is formed of a superconducting loop with three Joseph- 
son junctions. Two of the three junctions are of equal 
size, with Josephson energy Ej and charging energy E c , 
while the third junction is smaller by a factor a. Fig- 
ure^, shows a representation of the energy levels versus 
the value of the external magnetic flux in the loop for 
a set of typical parameters. The qubit quantum state 
can be represented as a superposition of two basis states 
that are persistent current states in the loop, with oppo- 
site values. Away from the symmetry point (fr q t, — <&o/2 
the energy eigenstates are almost equal to the persistent 
current states. When Q q b approaches $o/2, the energy 
eigenstates are superpositions of the basis states and for 
& q i, = $0/2 they become the symmetric and antisymmet- 
ric combinations of the basis persistent current states. A 
representation of the expectation value of the current in 
each energy eigenstate is given in Fig. QJd. 

The DC-SQUID is characterized by the gauge- 
invariant phase variables across the two Josephson junc- 
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tions, denoted by 71 and 72. The two variables are 
connected through the fluxoid quantization condition, 
71—72 = ^ 27r$ sg /$0i where $0 is the flux quantum 
and Q sq is the total flux in the SQUID loop. The flux 
in the SQUID loop contains an external component & x 
and a self generated component, which can be neglected 
for the typical parameters we will discuss. With this as- 
sumption, 71 — 72 = —2Ttf x , where f x — <fr x /<&o and the 
SQUID can be described as a single Josephson junction 
with a critical current given by I c (f x ) = 2J c q| cos (irf x ) | 
for a symmetric SQUID (I c i = I C 2 = I c o)- For symmetric 
qubit-SQUID coupling (when the two SQUID branches 
have mutual inductances to the qubit loop with opposite 
value) the flux generated by the SQUID in the qubit loop 
is given by M s I C i rc , where M s is the inductance between 
the qubit and the SQUID loops and I C i rc = (I\ — I2) /2 
is the circulating current in the SQUID loop. The latter 
is given by 



I arc = IcO Sin (tT f x ) COS"/ e , 



(1) 



where j e = (71 + 72) /2. The current and the voltage of 
the SQUID are related to the variable -f e through the two 
Josephson relations: 



and 



I = Ic {fx) sin j e , 



= £0 rfTe 

2tt dt 



(2) 



(3) 



From Eq. [21 and it follows that in the linear approx- 
imation (sin (7 e ) — 7 e ) the SQUID behaves as a linear 
inductor with the Josephson inductance 



2nl r 



(4) 



If an AC current is injected in the SQUID at frequency 
uq with a small amplitude J ac , the voltage across the 
SQUID has the amplitude V = 2nvQL,jI ac = $QVoI ac / I c . 
The maximum voltage across the SQUID is very small 
(~ 2/iV for vq = 1 GHz); very low noise amplification 
is necessary to detect such a voltage in a short time. 
Increasing the value of vq will result in a proportional 
increase in the value of the maximum AC voltage. How- 
ever, from Eq. ^and [21 it follows that when the SQUID 
current varies at frequency vq, the circulating current 
contains a significant frequency component at 2vq and 
additional higher harmonics for strong driving in the non- 
linear regime. The flux generated by this circulating cur- 
rent in the qubit loop can cause transitions between the 
qubit energy levels if the harmonics of the driving cur- 
rent are close to the qubit energy levels splitting. With 
typical level splitting of 1 — 20 GHz, the value of vq is 
limited to < 1 GHz. 

The relative change in Josephson inductance when the 
qubit evolves from the ground to the excited state is given 
by 



SLj 




SI C 


L.J 




Ic 



27 r|tan(7r/ x )|^, 



where it was assumed that the measurement is performed 
at a bias flux in the qubit away from <E>o/2, so that the 
expectation value of the qubit current in each energy 
eigenstate approaches in absolute value I p - the maxi- 
mum value of the persistent current (see Fig. QJ?) . The 
typical values for M s and I p limit the value of SLj/Lj 
to a few percent. If the SQUID is driven with a con- 
stant AC current, the maximum difference in AC voltage 
corresponding to a qubit state change, from the ground 
to the excited state, is ~ ^ v SLj/Lj. This can be in- 
creased if the SQUID is placed in a resonant circuit and 
the frequency uq is taken close to the circuit resonance 
frequency (see Fig. [TJ;). A limit on the quality factor Q 
will be set by the fact that the response time of the res- 
onator, Q/u)q, has to be smaller than the intrinsic qubit 
relaxation time, which is in the microseconds rang o 7 i 16 . 
When Q > 2Lj/6Lj the two circuit resonance peaks, 
corresponding to the different qubit states, are separated 
and a further increase of Q will not contribute to an in- 
crease in the AC voltage difference. 

The above considerations show that, given the typical 
qubit energy level splitting and relaxation time, the con- 
straints on the circuit parameters are vq < 1GHz and 
Q < 100. 



III. DETECTOR RESPONSE FUNCTION 

In this section we analyze the DC-SQUID in the in- 
ductive mode as a flux detector. We consider the case of 
moderate AC driving, when the SQUID behaves as a lin- 
ear inductor. The function describing the conversion of 
flux in the SQUID loop to AC voltage is determined for a 
general type of circuit in which the SQUID is embedded. 

If the magnetic flux in the SQUID loop varies in time, 
the relation between the transport current and the volt- 
age across the terminals of the SQUID is given by 



I(t) = 



1 



LAt) 



V{t')dt', 



(6) 



where Lj (t) is the time dependent Josephson inductance. 
Let us consider 



(7) 



(5) 



where a (t) parameterizes the variations of the magnetic 
flux in the SQUID. The time dependent Josephson in- 
ductance Lj(t) can be represented as the parallel combi- 
nation of the inductances Ljq and Ljo/a(t) (see Fig. 01). 
In the case of the qubit measurement, a(t) describes the 
dynamics of the qubit generated flux. The extreme values 
of a(t) in this case correspond to the qubit in a clockwise 
or anticlockwise persistent current state and are given 
approximately by ±7rtan(7r/ a; )M/p/$o (see Eq. [5j . We 
assume \a(t) | << 1, consistent with the usual value of 
the flux generated by the coupled qubit in the SQUID 
loop which is of the order of 1% of $ 2il2il&. From Eq.[f)] 
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(a) 



l c (t):c«— ►L J (t). 



(b) 




L J0 /«(t) 




- I — I - 



! L_in ! 




L J0 /«(t) 

FIG. 2: (a) The time dependent Josephson inductance can be represented as a series combination of the inductors Ljo and 
Ljo/a(t). (b) Schematic representation of the circuit in which the DC-SQUID is inserted (see the text for explanations). 



and[7| it follows that the current in the SQUID can be 
written as 



I = I +h, 



with 



and 



Io(t) 



h(t) = 



-JO 



«(*) 



■i JO 



V(t')dt' 



V{t')dt'. 



(8) 



(9) 



(10) 



As discussed in section [n] the measurement of the 
Josephson inductance is more efficient if the DC-SQUID 
is integrated in a resonant circuit. In this case the circuit 
is driven with an AC current source at a frequency close 
to the resonance frequency and the AC voltage is suitably 
amplified. The output AC voltage depends on a(t). Since 
there is a certain freedom in the design of the resonant 
circuit, we calculate here the dependence of the output 
voltage on a(t) for a general type of circuit in which the 
DC-SQUID is embedded. We consider a linear network 
with three ports that contains Ljo (the constant com- 
ponent of the SQUID inductance), the impedance of the 
driving current source, the amplifier input impedance, 
complemented by other linear circuit elements. The port 
Pi has its terminals across the current source. The port 
P m is connected across the Josephson inductance Ljq. 
Finally, the port P a has its connections at the input of the 
amplifier. Three elements are connected to the ports Pi, 
P a and P m , respectively: the ideal current source Ii(t), 
an ideal voltage amplifier and the inductance Ljo/a(t). 
A specific electrical circuit described in the way indicated 
here is shown in Fig. 

The relation between the current and the voltage at 
port P m is determined by the inductance Ljo/a(t): 



-Im(t) 



a(t) 



J JO 



V m (t')dt' 



(11) 



which is Eq. 1101 with changed sign in order to preserve 
the sign convention for the 3-ports network. The voltage 



at the port P a , a — o,m can be written as 
V a (t) -- 



Z ai (t')Ii(t-t')dt'+ / Z am (t')I m (t-t')dt' . 
'o Jo 

(12) 

Here Z a p{t) is the impedance matrix for the 3 port net- 
work, with a, (3 — i,m,o. We assume that the AC 
driving current is h(t) = I e cos(wot) and the flux vari- 
ations in the DC-SQUID loop (see Eq. [7J are described 
by a(t) — Re(cto e~xp(—iujt)) . Equation 1111 implies that 
V m and I m have components at frequencies wo + nw, with 
n being an integer. The voltage V m can be written as: 

Vm(t) = ^(^ e -»(-o+^)t + Y^hH e iWo+nu,)ty ^3) 



A similar expression for I m can be written if V is replaced 
by / in expression 1 131 

From Eq. ll2l written for a — to, one obtains the Fourier 
components of V m (t) as a function of the Fourier coeffi- 
cients of I m (t) (note that Ii(t) is imposed and has fre- 
quency components at ±loq). These values can be re- 
placed in 1111 and the terms corresponding to the fre- 
quencies loq + nuj are separated. Using the equations 
corresponding to n = 0, 1, —1 in expansions of the form 
1131 and neglecting the terms I m ,2 and I m ,-2 , the values 
for l m ,0i Im,i and I m ,-i can be obtained in lowest order 
in |oo|. Using these values in Eq. ^Jfor (3 = o leads to 
the following expression for the components of the output 
voltage at frequencies loq + u> and ujq — u: 



and 



ao 



ILOqL 



.10 



i{LUo)^-Z om (uj + uj) (14) 



-Z r . 



II) 



i( UJ o)-^Z 



(ujq-uj) (15) 



The expressions El arL d El f° r the up- and down- 
converted voltage at the output of the circuit are pro- 
portional to the amplitude of the driving current I e and 
to the amplitude of the flux modulation |ao|. These ex- 
pressions are usable only when the driving conditions (I e 
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and uo) are such that the maximum amplitude of the 
current in the SQUID is not close to the SQUID criti- 
cal current. Besides the up- and down- converted com- 
ponents V ,i and V ,-i, the output voltage contains a 
strong component U 0j o at frequency u) Q . V depends only 
quadratically on |ao| and thus it cannot be used for an 
efficient detection of the flux. 



The statistical properties of (t) are thus determined 
by 7^ (i) . From ^| it follows that for the case when the 
phase oscillations amplitude is large compared to the typ- 
ical values of 7 en (i), the most important contribution 
will be the mixing term "f e ,coh{t)"fe,n(t) ■ This results in 
frequency conversion of the circuit noise. 

The analog of equation^] for the quantized system is 



IV. QUBIT-SQUID INTERACTION 



¥e(t) = le,coh(t) +%n{t) 



(20) 



In a basis formed of two persistent current states, the 
Hamiltonian of the flux qubit can be written asii 



"■qb — T, a z + TTCx 



(16) 



where &i, i — x,y,z, have the Pauli matrices repre- 
sentation. The coefficient of the first term in ^| is 
e = 2I p - $ /2), where $ gt is the flux in the 
qubit loop. The maximum persistent current I p and 
the minimum energy level splitting A (see Fig. and 
b) are parameters fixed by the qubit junctions design. 
The average flux induced in the SQUID loop by the 
qubit is MI p (a z ). The flux dependent term in the 
energy of a SQUID is the Josephson energy given by 
-2$ /(27r)/ c0 cos(7 e )cos(7r<l) ;! . (Z /(f>o). The total flux 3> sg 
in the SQUID loop contains the external flux & x and 
qubit induced flux MI p (a z ). It follows that the interac- 
tion Hamiltonian can be written as 



H c = MIpIco cos(%) sin(irf x )a z 



(17) 



where we assumed that the flux generated by the qubit is 
small and thus a linear approximation could be used. A 
rigorous derivation of the interaction term in the Hamil- 
tonian for a coupled DC-SQUID and a three Josephson 
junctions qubit, assuming a SQUID with a small self in- 
ductance and using the two level approximation for the 3 
Josephson junctions loop leads to the same result as 1171 
If OH is compared to QU with e/2 = I p (<P qb - $ /2), 
it becomes clear that the back-action due to the mea- 
surement is described by an equivalent flux operator, ex- 
pressed as = MI C i rc . I c i rc is the operator correspond- 
ing to the circulating current in the DC-SQUID and is 
given by I circ = 7/cos(7 e ), where 2/ = I ca sm(Trf x ) (see 
Eq.P). Classically, qubit decoherence can be understood 
as a result of the fluctuations in the flux bias, due to the 
SQUID. If 7 e is treated as a classical variable, its time 
evolution is given by 



7e(*) = le,coh(t) + 7 e ,„(t) 



(18) 



where 7 e , co h(t) = Re (7 e0 exp (-iuj t)) with 7 e0 = 
I e Zi m (loq) I (iuJo) the response to circuit driving and 
le,n{t) is a random term, corresponding to e.g. ther- 
mal fluctuations. The "classical" flux is given by $>b(t) = 
MIf cos(7 e (t)), which can be approximated by 



S t (t) = MI f 1 - 



llit) 



(19) 



in which 7e(t) is the phase operator in the interaction 
representation with respect to the qubit-SQUID interac- 
tion, which is thus equivalent to the Heisenberg repre- 
sentation for the SQUID system. ^ e ,n{t) is the phase 
operator representing the intrinsic evolution, in the ab- 
sence of circuit driving. The first term on the right hand 
side of Eq. [201 corresponds to the evolution due to driv- 
ing, and it is the same as the first term on the right hand 
side of Eq.HHI 

The essential feature of the interaction Hamiltonian 
given by 1171 is that the coupling to the external phase 
operator does not have a linear part. Recent work on the 
influence of nonlinear coupling of the noise on the evo- 
lution of a two level system has been done22 motivated 
by results reported by Vion et al&, where long coherence 
times were obtained for the operation of a qubit at set- 
tings where the energy level separation was insensitive 
in the first order to external noise. In the second order 
approximation, the back-action flux noise is described by 



<f b = MI f [ 1 - f 



which can be separated in three parts as: 



(21) 



61 



MI 



1 - 



7 2 

/ en 



(t) 



"1e,coh 



2/2 

(22) 

The first term on the right hand side of!22lcan cause tran- 
sitions between the qubit energy levels. As we discussed 
in Sec [Hi resonant transitions occur when the qubit en- 
ergy levels splitting is close to the harmonics of the driv- 
ing AC frequency, and in particular to the second har- 
monic. The time average of the first term in 1221 depen- 
dent on the amplitude of the AC driving current, will be 
considered a part of the qubit flux bias The effects of 
the second term in [23 were analyzed by Makhlin et al22k 
for ohmic and 1/lu type spectral densities. In this paper 
we focus on the calculation of decoherence determined by 
the third term in [22 In the second order perturbation 
theory, used for the calculation of the decoherence rates 
in the next section, the contributions from the different 
terms in 1221 can be treated independently. 
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V. CALCULATION OF THE DECOHERENCE 
RATES 

In this section we calculate the relaxation and dephas- 
ing rates of a flux qubit during the measurement by a 
DC-SQUID in the inductive mode. It is assumed that the 
external qubit flux <$> q b is fixed. However, the calculations 
can be extended to include the case of the measurement 
performed during induced Rabi oscillations or other con- 
trol sequences^S*. The model Hamiltonian used for the 
combined system qubit-SQUID is: 



H(t) — H q b + H c + HsQuiD{t), 



(23) 



where H q (, is the qubit Hamiltonian given by 1161 H c is 
the interaction term given bv !17l and HsQUiDit) is the 
SQUID Hamiltonian, which is time dependent due to 
driving. If a transformation is made to the qubit energy 
eigenstates, the first two terms in 1231 become: 



Hqb = 



A 2 



and 



H c = MI p If cos(%) {cos(6)t z - sin(9)f x ) 



(24) 



(25) 



where the fj, i — x,y, z are Pauli matrices in the energy 
eigenstate basis and tan(0) = A/e. In the interaction 
representation with respect to the qubit-SQUID interac- 
tion, the operators r\ evolve in time according to 

T l I (t)=A ij (t)r j (26) 

with the matrix A given by: 

cos(tJoii) — sin(u;oii) 
A(i) = ( sin(woit) cos(w ii) ] (27) 
1 



evolution of the operators rf- in the Heisenberg picture 
is obtained using time-dependent second order perturba- 
tion theory. As explained in Section llVl the relevant part 
of the interaction Hamiltonian (Eq. I25fl for the calcula- 
tion of decoherence is given in the interaction picture by: 



Hl(t) = -MI p I neiCoh {t)%, n {t){cos{6)Tl - sin(0)fi). 

(28) 

The evolution of the operators f t - , which allows to de- 
scribe the qubit operators expectation values if the initial 
state is known, is given by: 



1 



t>ti>t 2 >0 



In the following, the second term on the right hand side 
of Eq. |2HI is calculated. For the initial calculation we 
assume the most general interaction Hamiltonian with 
linear coupling to the bath, which can be written as 



Hl{t)= //(*)*/(*)> 



(30) 



i—x,y,z 



where /j are bath operators (note that the interaction 
representation is used in !3UI) . In the end the form of the 
operators //(t) corresponding to our case, as given by 
Eq. will be considered: 



flit) 
flit) 



MIpIf sm(9)j e , coh (t)^ en (t) 


-MI p If COs(0)j e ^ oh (t)% n (t). 



(31) 



in which u)qi = -\/e 2 + A 2 /fi is the frequency correspond- 
ing to the qubit energy level separation Ve 2 + A 2 . The 



If the commutation and anti-commutation relations for 
the fj operators are used, Eq. [25] results in 



f(*) = 



Aij(t) 



^ J J dhdt 2 (d+itut^n - O+ih^fj - e jlm 6i m {t x M) 



(32) 



t>*l,*2>0 



In the last expression: and 

6£(ti,ta)= J2 Ai fe (-t 1 )C±(t 1 ,fe)Ay(t 2 ) (33) C^h) = ^(/ fe J (ii)//(i 2 ) - fl^fUh)). (35) 

k,l=x,y,z ^ 

with The last two expressions are symmetrized and anti- 

„ i A j -j -j -j symmetrized products of operators at different times. 

C+(ii,t 2 ) = ^(f k (ti)fi {t 2 )+fi (t2)f k (ti)) (34) 

Their expectation values calculated for a thermal equilib- 
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rium state are connected with the linear response func- 
tions by the fluctuation-dissipation theorem 2 ^. Note 
that, to obtain |221 the integral in [2HI was extended to 
the region t\ < t 2 because the integrand is symmetric 
under the interchange of t\ and £2- 

We assume that initially the qubit and SQUID states 
were separable and the SQUID is described by the ther- 
mal equilibrium density matrix. For the case of the cou- 
pling Hamiltonian given in Eq. 1281 the relevant correla- 
tion functions are: 



ci(t 1 ,t 2 ) = (-{^(t 1 ),^(t 2 )} + ) ( 



and 



£7. ^2) 



(^{¥ e (ti),¥ e (h)}-)o, 



(36) 



(37) 



where +/- denote the anticommutator / commutator, and 
the expectation value is taken for the SQUID thermal 
equilibrium density matrix. From l3*21 - l37l we see that the 
time evolution of the operators fi, i = x, y, z, depends on 
their expectation values for the initial qubit state and on 
a 2 dimensional integral involving the expectation values 
of operators of type and |33 

The interaction between the qubit and the measure- 
ment DC-SQUID has the consequence that the qubit 
quantum state becomes a mixed state. In general one 
distinguishes between energy relaxation, corresponding 
to a change in the qubit energy expectation value, and 
dephasing, corresponding to randomization of the phase 
of a coherent superposition of energy eigenstates^. To 
calculate the energy relaxation, we determine the tran- 
sition rates between energy eigenstates by determining 
the evolution of (f~(t)). If in this calculation the ini- 
tial qubit state is chosen to be the ground or the ex- 
cited state, these rates will represent the absorption and 
emission rates respectively. To calculate the dephasing 
rate, we determine the decay of the expectation values 
(l/2(f^(i) ± iTyit))), with the qubit initial state being 
a t x eigenstate. 

The calculation of the integral on the right hand side 
of 13 21 involves a product of the functions cos(woiii,2) and 
sin(£Joiii,2)> resulting from the expression for the free 
evolution matrix A (see Eq. 1271 and of the functions 
cos(wo^i,2) and sin(woii,2), resulting from the time de- 
pendence of the coupling operators //(t) (seel5T)l. The 



correlation functions appearing in 1361 and 1371 only de- 
pend on the time difference t\ — t^. For times that 
are long compared to 27r /loq and 27r/woi the relax- 
ation (decay of (t^ (i))) and the dephasing (decay of 
(1 1 '2{t^ {t)±iTy it))) ) can be described as an integral of 
the product of the Fourier transform of one of the spec- 
tral functions C^itx,^) and a weight function that has 
a width depending on the integration time t. This weight 
function is given by: 



dt\dt<xe 



-iu(ti-t2) 



(38) 



0<*i,i 2 <t 

and has the property lim W(u>,t) 



5(lo) . The expres- 

sions below for the relaxation and dephasing rates are 
given assuming that the spectral density of the circuit 
noise, given by the Fourier transform of the correlation 
functions 1361 and 1371 does not have significant variations 
over the frequency range where the weight function has 
substantial values. In this case, the relaxation and the 
dephasing of the qubit state are proportional to the time 
from the beginning of the measurement. 

The transitions rates between the two energy eigen- 
states depend on the initial state. The transition rate 
from the excited state to the ground state T± (emission) 
and the transition rate form the ground state to the ex- 
cited state T| (absorption) are given by: 

Tx = ^ sin 2 (6»)fc 2 (7 e o)[5 7e (woi + w ) + 5*7.(^01 - ^0)] 

(39) 



and 



(40) 

in which relations 

fc(7 e o)=MJ p 7 / | 7e o| (41) 

is a measurement coupling factor and 

S 7 >)=S+(w)-i£-(w). (42) 

S^ e are the Fourier transforms of the correlation func- 
tions (i, 0) given by an d U%7\ The corresponding 
expression for the dephasing rate is 



T = ^sin 2 (0)fc 2 ( 7e o)[S+(woi +^ ) + 5-(w i ~ wo)] + ^cos 2 (0)fc 2 ( 7eO )S+M. 



(43) 



r 



VI. DISCUSSION 

In this section we discuss the results of the calcula- 
tions of the parameters characterizing qubit decoherence 



and we analyze a practical circuit which can be used for 
single-shot readout of a flux qubit. 

We start with a discussion on the emission and ab- 
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sorption rates, along the lines of similar analysis done 
for charge^ and charge-phase^ qubits. Both and Tj 
are proportional to sin 2 (0), due to the fact that the op- 
erator t x (see I25|l causes transitions between the energy 
eigenstates. The difference between the two rates given 
by [2HI and 0U| is due to the last term in the integrand in 
expression 1321 connected with the fact that a commuta- 
tor is non- vanishing, so it can be attributed to quantum 
noise. Let p g and p e be the probabilities for the qubit 
to be respectively in the ground and in the excited state. 
The time evolution of these probabilities is determined 
by the rate equations: 



d Pg 
dt 

dp e 
dt 



r. 



-Pa 



'-Pa 



r 



-Pe 



(44) 



and the normalization condition p g +p e — 1. Since and 
Tj in 1391 and 1401 describe the decay of (f 2 ), they appear 
divided by 2 in EH The polarization P(t) = p g {t) — p e {t) 
tends to the equilibrium value 



with a relaxation rate 

r,. 



r i+ r T 



(45) 



(46) 



The spectral densities of the symmetrized and anti- 
symmetrized correlation functions 1361 and 1371 depend on 
the impedance at the port P m (see Section 1111}) as: 



and 



• coth 



-iS (w) 



2k B TJ ' 



Rk 



8tt Re(Z mm (uj)) 
oj R K 



(47) 



(48) 



where Rk = h/e 2 (see DevoretSl). Given the relations 
EHHHIEII|andE2]and the properties of S+ (u) and S~ (uj) 
to be respectively even and odd functions, the relaxation 
rate can be written as 

T, = ^ sin 2 (0)fc 2 ( 7eO )[5+ (woi + lu ) + S+ (woi - wo)]. 

(49) 

For a flux qubit coupled to a DC-SQUID biased with a 
constant current van der Wal et ali* found that qubit 
relaxation is proportional to (o>oi). Our results show 
that because of driving the SQUID with an AC bias cur- 
rent at frequency luq the qubit relaxation rate is propor- 
tional to the sum of (cJqi +^o) and 5! + (uj qi —uiq), and 
multiplied by the coupling factor of Eq. jjfl] In practice 
c^q << u>oi, which implies that for a spectral density of 
the noise which is reasonably flat at large frequencies we 
can take 5+ (w i + ^o) ~ >S+,(woi - ^o) ~ 5+ (o>oi) and 



our results is not significantly different of the case of a 
DC current biased SQUIDil 

From E3 E3 and El it follows that 



^7e(-^) 



-phu 



(50) 



where = l/(k B T). For a qubit coupled to a DC-SQUID 
biased with a constant current, which is the case analyzed 
by van der Wal et a/:—, r j and T| are proportional to the 
spectral density of the noise at the frequencies u>oi and 
—woi , respectively. For that case Eq. 1501 implies that 
T T /rj. = e^" 01 . This is the detailed balance condition 
and implies that in a stationary situation the qubit is 
in thermal equilibrium with the environment at temper- 
ature T (see also the results of Schoelkopf et alm§>). In 
contrast, the relations 1391 and I4TH show that for the case 
analyzed here the detailed balance condition is in general 
not satisfied. 

The dephasing rate in 1431 can be written, if 1461 is used, 

as 



L 4> = : \- L d» 



where the pure dephasing rate is given by 



r; = -cos 2 (#)fc 2 ( 7e0 )s+M. 



(51) 



(52) 



The factor cos 2 (8) is due to the coefficient of the operator 
f 2 in 123 Dephasing is a result of random modulation of 
energy level separation due to noise in the SQUID circu- 
lating current. The fact that the SQUID is driven with an 
AC current has the consequence that the pure dephasing 
rate depends on noise at loq, which is qualitatively differ- 
ent of the result obtained by van der Wal et alH. For the 
radio-frequency Bloch-transistor electromete»A& a similar 
contribution of the converted noise to back-action was 

We compare the pure dephasing rate FJ given bv 1521 
with a similar contribution due to the second term in Eq. 
I21las calculated by Shnirman et alM, that we denote by 
ri. We consider the simple case where the DC-SQUID 
is shunted by a resistor R s h , corresponding to ohmic dis- 
sipation, when the result of Shnirman et al. can be used. 
The following relation is valid for ujqLjq << R s h'- 



p* 



coth 



^0 \ 7e0 RshPK 



ksT J 2 {loqL 



JO J 



(53) 



For the case loqLjq « Rk and fkvo ~ fcsT, the dephas- 
ing rate is dominant even at small SQUID driving 
amplitudes. 

The reliable measurement of the qubit state requires 
that the AC voltage at the output of the circuit is aver- 
aged for a long enough time, such that the noise due to 
the amplifier is less that the difference between the volt- 
age values corresponding to the two qubit flux states. We 
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define the discrimination time as the time necessary to 
have a measurement signal to noise ratio equal to 1. It 
is thus given by 

Sy (cjp) 

dlscr ~ (AV qb ) 2 ' [ ' 

where Sy(u>o) is the spectral density of the voltage noise 
and AV q b is the difference in the output voltage values 
corresponding to the two qubit states. The value of AV q b 
is proportional to 7 e o- The discrimination time Xdj scr , 
the relaxation time T r — \/T r and the dephasing time 
Ttj, — 1/r^, are inversely proportional to 7^0 (see 1461 with 
|2U and and I43|l . Increasing the amplitude of the AC 
driving leads to a decrease in the discrimination time. 
However, this is accompanied by a proportional decrease 
of the qubit decoherence times T r and T^. This illus- 
trates the tradeoff between obtaining information about 
a quantum system and the state disturbance. The mea- 
surement is efficient if the ratio T r /Tdi scr is large. The 
ratio T r /Tdiscr does not depend on the amplitude of the 
AC driving. However, a fast measurement is necessary if 
we take into account the fact that, besides the measure- 
ment back-action, there are also other sources of deco- 
herence that will increase the total relaxation rate. 

We now analyze the measurement of a flux qubit using 
our particular SQUID embedding network presented in 
Fig. 01. The bias resistor Rb has the purpose to increase 
the impedance of the current source. The inductor L s is 
a small stray contribution, unavoidable in the design of 
the circuit. The combination of the capacitors C\ and 
C2 is an impedance transformer that will increase the ef- 
fective impedance of the amplifier input, at the cost of 
a division of the total voltage across the inductors; they 
also provide the capacitive part necessary to create a res- 
onant circuit. The DC-SQUID has Josephson junctions 
with a critical current I c o = 200 nA. The external mag- 
netic flux in the SQUID loop corresponds to f x — 3.35, 
resulting in a critical current I c = 187 nA. The measured 
persistent current qubit has I p = 300 nA and A = 5.5 
GHz. Figure^, and b show plots of the energy eigenval- 
ues and persistent current expectation value versus bias 
flux for these qubit parameters. If the mutual induc- 
tance between the qubit and the DC-SQUID is M = 40 
pH, the relative change in Josephson inductance is given 
by a = 3.4 %. A plot of the expression S Jb (ui) given 
by 021 with El an d EH is shown in Fig. |3t>, assuming a 
temperature T = 30 mK. 

To calculate the discrimination time given by Eq. 1541 
we assume that the voltage noise is dominated by the 
voltage amplifier. We assume that a low noise cryo- 
genic amplifier with a noise temperature of 4 K is used^S. 
Equations 1141 and 1151 allow the calculation of AV^. We 
assume that to ~ 0, since qubit relaxation is slow com- 
pared to the detector bandwidth (which will be confirmed 
by our calculation of the relaxation time) and we choose 
the value loq = 2-7T x 672 MHz that gives a maximum am- 
plitude AV q b = 189 nV for an AC driving current such 
that the amplitude of the SQUID phase oscillations is 




0.490 0.495 0.500 0.505 0.510 



FIG. 3: (a) Schematic representation of the measurement 
circuit, with notations according to Sec. 11111 The values of 
the circuit elements are Rb = 4.7 kfl, Z s = Z a — 50 Q, Ljo = 
1.76 nH, L s = 0.18 nH, C x = 60.7 pF and C 2 = 60.6 pF. 
(b) Representation of the Fourier transform of the correlation 
function for the SQUID phase operator versus frequency, (c) 
Plot of the measurement discrimination time (continuous line) 
and qubit relaxation time T r (dashed line) versus qubit bias 
flux. The measurement is performed with an amplitude of 
the AC driving such that 7 e o = 0.5 at a frequency vo — 672 
MHz. 



7 e o = 0.5. The discrimination time is plotted in Fig. 
The discrimination time increases when the qubit bias 
flux $ 9 approaches $o/2, because the difference between 
the expectation values of the qubit current for the two 
energy eigenstates decreases (see Fig. J2Jd). 

The relaxation time is calculated using the expres- 
sion 021 an d plotted in Fig. ^Ja for the chosen operating 
frequency ujq. The relaxation time away from the sym- 
metry point $ 9 = $0/2 increases as a result of both the 
decrease in the transverse coupling term sin 2 (#) and the 
decrease in the real part of the impedance Z mm away 
from the resonance peak (see Eq. 0^1 and I47|) . Over a 
wide range of parameters the relaxation time is consider- 
ably higher than the discrimination time, which allows for 
very efficient readout of the qubit state. Using a SQUID 
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amplifier— with a noise temperature less than 100 mK 
would allow for the reducing the discrimination time by 
more than one order of magnitude. 

The measurement of the qubit state can be performed 
by applying the AC current to the SQUID for a time T m 
and by measuring the average AC voltage during this 
time interval. Note that the readout does not have to be 
performed at the same qubit bias flux where qubit manip- 
ulation prior to measurement is performed. It is possible 
to perform operations on the qubit at <S> q i, = $o/2, where 
the qubit is insensitive to magnetic flux fluctuations. Af- 
terwards the flux in the qubit can be changed adiabat- 
ically to a value where the two energy eigenstates have 
sufficiently different values of the persistent current to al- 
low discrimination by the measurement. If the measure- 
ment is performed at $ g f, — 0.497 $o where Tdiscr = 150 
ns and T r = 50 /is, a measurement time T m = 500 ns 
ensures a measurement fidelity larger than 80 % and the 
qubit relaxation is negligible during this time interval. 

The dephasing time depends on the Fourier transform 
of the symmetrized correlation function at the frequency 
of the AC driving. It follows from 03 that 5+(wo) is 
large, because u>q has to be close to the resonance fre- 
quency of the circuit for efficient state readout. Even a 
small amplitude of the AC signal can cause significant de- 
phasing of the qubit. During qubit manipulation, when 
no measurement is performed, the SQUID AC driving 
current has to be suppressed very strongly. For operation 
at a qubit energy level splitting cjqi = 2A the decoher- 
ence time due to the SQUID is 10 fj,s if the amplitude of 
the phase oscillations is |7 e o| = 0.003. 



The continuous nature of the flux detection makes this 
readout method suitable for fundamental studies of the 
dynamics of the measurement process. Further analy- 
sis will be necessary for understanding the dynamics of 
the coupled qubit- SQUID system and for an evaluation 
of possible direct observation of qubit coherent evolu- 
tion, similar to the situation described by Korotkov and 
Averin— . 



VII. CONCLUSIONS 

In this paper we analyzed the DC-SQUID in the induc- 
tive mode as a readout method for superconducting flux 
qubits. We characterized the response function of the 
DC-SQUID as a flux detector. We described the back- 
action of the measurement circuit on the qubit. The re- 
laxation and dephasing rates are proportional to circuit 
noise at frequencies that are shifted by the SQUID AC 
driving frequency, which is a result qualitatively differ- 
ent of the case of a the measurement done with a switch- 
ing DC- SQUID. For a realistic measurement circuit, we 
found that single shot measurement of a flux qubit is 
possible. 
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